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Abstract 

In the solar atmosphere, the chromospheric and coronal plasmas are much hotter 
than the visible photosphere. The heating of the solar atmosphere, including the 
partially ionized chromosphere and corona, remains largely unknown. In this paper we 
demonstrate that the ions can be substantially heated by Alfven waves with very low 
frequencies in partially ionized low beta plasmas. This differs from other Alfven wave 
related heating mechanisms such as ion-neutral collisional damping of Alfven waves 
and heating described by previous work on resonant Alfven wave heating. In this 
paper, we find that the non-resonant Alfven wave heating is less efficient in partially 
ionized plasmas than when there are no ion- neutral collisions, and the heating efficiency 
depends on the ratio of the ion-neutral collision frequency to the ion gyrofrequency. 
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I. INTRODUCTION 



The chromosphere, the region between the solar surface and the corona, is permeated by 
low-frequency (0.0014 Hz ~ 0.0079 Hz) Alfven waves with strong amplitudes (2.6 kilometers 
per second in amplitude) [I| while the plasma beta value(/3 = (v p /va) 2 , where v p and v& are 
the thermal speed and the Alfven speed, respectively) in this region is low (where /3 reaches 
a minimum with /3 < 3 x 10 ~ 2 ) [2J,[3] - The chromospheric and coronal plasmas are much 
hotter than the visible photosphere. The heating mechanisms in these regions, however, 
have not yet been fully understood [HEJII] . Alfven waves have long been considered to play 
a crucial role in heating of plasma in these two regions and in magnetic fusion devices [5J|B]. 
Numerous theoretical and experimental papers have been published to investigate resonant 
heating of ions by Alfven waves [THTTj. In these works, the cyclotron resonant condition 
is necessary for ion heating by the Alfven waves, and in general the frequencies of the ap- 
plied Alfven waves are comparable to the cyclotron frequency. However, the heating of ions 
by low-frequency Alfven waves and related problems have triggered great interest in recent 
years [T2"H2"2"] . In these works, the ions can be heated by Alfven waves via nonresonant 
interactions, which means u <C Q q , where u is the frequency of the Alfven waves and Q q 
denotes the gyrofrequency for ion species labeled q. Researchers have used various simula- 
tion and theoretical methods in order to validate the heating mechanism of low-frequency 
Alfven wave interactions; these include test particle approaches [121 dHJ EJUS! ED] , two-fluid 
magnetohydrodynamic (MHD) model (nonlinear parametric interaction of Alfven waves with 
slow/fast waves) [H], hybrid simulations [I2JI22], an d kinetic theory [T8||2"2"]. Importantly, 
several of these works demonstrate quantitatively similar results even while employing dif- 
ferent approaches. For example, the numerical and analytical results adopting the heating 
mechanism proposed by Wang and Wu (discussed later in this paper) [T71HB] are quantita- 
tively in accordance with the work of Voitenko and Goossens who examine cross-field heating 
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by low-frequency kinetic Alfven waves (KAWs) [T6] . 

Some recent works [231425] indicate that Alfven waves can propagate through the par- 
tially ionized solar chromosphere. They find that Alfven waves with frequencies below 0.6 
Hz will not be completely damped. More importantly, these works all show that Alfven 
waves with frequencies below 0.01 Hz are undamped by ion-neutral collisions in the solar 
chromosphere, and are therefore available for damping by other mechanisms, such as the 
mechanism described in this paper. 

In this paper, we demonstrate that the low-frequency Alfven waves propagating along 
the background magnetic field B = B i z , can heat ions even in partially ionized plasmas. 
It is important to note the heating mechanism in this paper is different from the previous 
works [2"3H2"6] . The heating process in the present work is due to a randomization of the 
spatial velocity distribution in the ion population, which is caused by the nonresonant Alfven 
wave interactions. We find that the heating process becomes less efficient than the situation 
with no ion-neutral collisions. Moreover, the nonresonant heating process is only effective 
for low-beta plasmas [T7] and Alfven wave frequencies lower than 0.6 Hz; the most efficient 
heating occurs when u < 0.01 Hz since low frequency Alfven waves are more abundant due 
to the damping of high frequency waves by ion- neutral collisions [2"3"H2"B] . 

II. ANALYTIC THEORY AND TEST PARTICLE CAL- 
CULATIONS 

We consider the Alfven waves have a spectrum and the dispersion relation can be described as 
u = kvA {va is the Alfven speed, u and k are the wave angular frequency and wave number, 
respectively). This relationship is still appropriate even when the plasma is partially ionized, 
as described in the following paragraphs. Without loss of generality, we consider left-hand 
circular polarization in this paper. The wave magnetic field vector ^B^, and electric field 
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vector 5~E W can be expressed as 

6B W = ^2B k (cos(j)kix - smcf) k iy), (2.1) 
k 

5E W = -^b x 5B W , b = 2° (2.2) 

c fio 

where i x and i y are unit directional vectors, <pk = k(v^t — z) 4- <^fc denotes the wave phase 
and (fk is the random phase for mode k. In the following, we pay attention to the protons 
only, whose equation of motion is described by 

= q% (SE W + ^ x (Bo + ^B^)) +miV in (u- v), ^ = v ( 2 - 3 ) 

where v is the ion velocity, u is the bulk velocity of a background neutral fluid, and v in is the 
frequency for elastic collisions between ions and neutrals. Ion collisions with electrons are 
neglected in Eq. fl2.3j) because m, 3> m e . The collision frequency responsible for momentum 
transfer between species i and n is defined as: 



m n 8k B T 

n n \ (Tin, (2.4) 



rrii + m n '"V vrm m 

with m in = (mim n ) / (rrii + m n ),o"j n the collisional cross-section for collisions between the 
two species. Here m; = m n because the plasma we studied in this paper is assumed to be 
entirely composed of hydrogen which leads to m in = mi/2. The wave phase speed v p h here 
is defined by [24] : 

v ph = - = v A > l-i (2.5) 

in which p n is the mass density of neutrals, ptot is the total mass density of the plasma, and 
v ni is neutral-ion collision frequency. Therefore, the relationship p n / Ptot < 1 is necessarily 
valid. As shown by De Pontieu and Haerendel, for waves with v < 1 Hz, the assumption 
w ^ v ni {u/vni ^ 1) holds throughout the chromosphere [23]. Furthermore, most theories 
for the generation of Alfven waves in the solar atmosphere predict typical frequencies below 
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1 Hz [21], which has also been observed pp. These allow us to simplify Eg. (12. 5ft to 



v ph = ~r = Va \ ~ v a ( 2 -6) 

Hence the relationship u = kvA is still appropriate. 

Defining v ± = v x + iv y , u±_ = u x + iu y , Vu = v z , Uu = u z and 5B U = J2k B k e~ lc ^ k ; we are 
left with 

^ + (ift„ + ^n)«l = i( V \\ - V A) ^2 fi fc e_l0fe + U iriU± (2.7) 

^ = - Im (v^Si k e irt>k ) + v in {u\\-v{), ^ = v\\ (2.8) 



where Qq = ^7 (the proton gyrofreguency) , Q k = Im{ ) denotes the imaginary 

part of its argument. As a first-order approximation, we can assume v\\ ~ wii(O), where 
f||(0) is the particles initial parallel velocity. The approximation is valid when ^ = 
is small enough and the freguencies of the Alfven wave are sufficiently low to ensure that 
I k (v\\(0) — va) I- For simplicity, we assume the bulk velocities of the cold neutrals 

u x (t) = u y (t) = u z (t) rs and ion- neutral collision freguency v in ps constant, which provides 
a lower limit for the amount of ion heating. Since the ions and background neutral fluid are 
both "cold" initially, it is reasonable to assume u(t) ~ due to the relatively short time 
scales (compared to the time reguired to heat the neutrals) considered in this paper and the 
high neutral fraction (detailed below). 

The ion gyrofreguency in the region of 800km ~ 1500km above the photosphere in solar 
chromosphere ranges from 1 x 10 5 ~ 5 x 10 5 Hz. The ion-neutral freguency in this region 
ranges from 5 x 10 3 ~ 5 x 10 5 Hz. These results are based on the solar atmospheric model 
VAL C [57] for the densities and temperatures, and a magnetic flux tube model with 1500G 
field strength in the photosphere and 10G magnetic field strength in the corona [2S]- We 
select two values of the ratio a (a = ^g-), 0.1 and 0.5, to represent the region around 
1000km above the photosphere. According to VAL C model [2Z] and the magnetic flux tube 



model [28] described above, when a = Ui n /Qo > 0.05, the ratio of neutral density to ion 
density A = n n /rii 1, which indicates that the amount of neutrals is much larger than 
that of ions. It is important to note that although the value of ratio A ^> 1 indicates that 
u in is large (with respect to Alfven wave frequency u; refer to Eq.f l2.4p ). the value of flo can 
also be large (with respect to u) due to the background magnetic field B (f2 = ^r), thus 
the value of ratio a can still be very small. 

In the following, we focus on the perpendicular velocity component due to the fact that 
the ion temperature increase is more prominent along the perpendicular direction than the 
parallel direction [TSltTTlIIH] ■ We acknowledge that the anisotropy between v± and v\\ should 
be alleviated after a sufficiently long period of time due to the collisional effects, however, 
this is not taken into account due to the relatively short time scales considered in this 
paper. Even over much longer timescales, where ions temperatures become isotropic, the 
fact that the ions are heated by Alfven waves remains unchanged. With the initial condition 
v± = vj_(0) and z = z(0), the solution of Eq. fl2.7p can be written as: 

v± = MOK^ ^"^ + UinU± f in - lQo) [1 - e-^e-"*] 

^0 + v in 

_ Zfc ^fc l v A ~V\\ (0)] (Aq + Mjn) r- ik ( VAt - z y itpk _ e i(fc z (o)- Vfe ) e -(iQo+i/i re )tj ( 2 .9) 

Here, we use the approximation that Q — k[vA — U|i(0)] ~ an d z = z(0) + vn(0)t. In 
order to verify that the analytical solution is correct, we set z^ n = 0. Then the Eqs. (12.31 
12.71 12. 8p are reduced to those found in previous work [T^ITOITH] where the analytical result 
is as follows: 

v± = v x (0)e- iQot - v A ^k B k e -ik{v A t-z)-i Vk + VA Lk B k e nkz{o)- Vk ) e -«iot (2.10) 

Bo B 

The relationships below between the kinetic temperature and velocity of protons are based 
on plasma which consists of an ensemble of protons, 

lJLmv 2 (t) J^rniVu(t) 

rxW-jE 1 ™ r «« = £^P (2 - n) 

i=l i=l 
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where T± is the perpendicular kinetic temperature, Tji is the parallel kinetic temperature 
and N is total number of the protons. 

We present the simulation results using test-particle calculations that build upon previous 
works [T31[T7] IT9~ | |2"0"]. The test-particle simulation will be valid in partially ionized plasmas 
when no ion-neutral collisional damping occurs. We discretize the Alfven wave number 
by kj = k min + (j - l) fcroa }:^" , for j = 1, J, where k min = k x = 1 X lO~ 8 Q /v A and 
kmax = kj = 5 x 10~ 8 Qq/va- This range of wave numbers implies that we are considering 
1 x 10~ 8 f2o < u < 5 x 10~ 8 Slo? so that the wave frequencies are much lower than the 
proton gyrofrequency and also will be low enough to guarantee that no ion-neutral collisional 
damping occurs [2"3"H2"5] . The field amplitudes of different wave modes are equal to each other, 
and they are constant. Here we set two values of SB^/B^ = ( ~ 0.05 and 0.12. The total 
number of test particles is 10 5 , which are randomly distributed during the time interval 
< Q()t < 2tt and along the spatial range < zQq/va < 3 x 10 9 . In this paper, we vary 
the ratio a (a = j^) between 0.1 and 0.5, which is well within the acceptable range for 
describing the solar chromosphere as determined by the previous works [2ZII2H]- The initial 
ion velocities are assumed to have a Maxwellian distribution with thermal speed v p , which 
is less than va to ensure that the cyclotron resonance condition cannot be satisfied. 

In Fig.l we present scatter plots in the v x — v y space which illustrate the process of 
particle heating. Here, the velocity is normalized to the initial thermal speed v p , which is 
set to be v p = 0.07v A {P = {v p /va) 2 = 4.9 x 10~ 3 1, which is well within the range of 
chromosphere beta value determined by Gary [3]). We compare nonresonant heating without 
and with neutral collisions in Fig.l. The results for the non-collisional case, where a = 0, 
shown in Fig.l (a), (b), (c), and (d), are in accordance with former results [17]. We find that 
ion-neutral collisions reduce the amount of heating over time and are responsible for creating 
a ring distribution. Since u(£) ~ is valid during the relatively short time scales considered 
in this paper as explained above, a balance is reached between energy gained via Alfven 
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Figure 1: Velocity scatter plots of the test particles in the v x — v y space in the case u in = 
(first line) and u in = 0.1Q (second line) at time fit = 0, fit = 7, fit = 20 and fit = 40, for 
input parameters (5B W 2 /B 2 ,v p /v A ) = (0.05, 0.07). 

wave interactions and lost from collisions with cold neutrals. Thus the velocity distribution 
evolves into a "thin" ring structure, with all of the ions evolving to a single speed. This 
result can also be analytically approximated, and will be discussed in detail in the following 
paragraphs. 

Fig. 2 shows the temporal evolution of the kinetic temperatures, where the results are 
based on two sets of input parameters (5B W 2 /B 2 , v p /va) = (0.05, 0.07) and (0.12, 0.07). The 
cases a = 0, 0.1 and 0.5 are represented by the dotted line, the solid line and the dashed 
line, respectively. The results again demonstrate that heating via nonresonant Alfven wave 
interactions is possible, and that even in the presence of neutrals heating still occurs and 



reaches a steady value. Fig.2 illustrates that the stronger Alfven wave amplitude 




will result in the larger amounts of heating, which is also indicated by Eq. fl2.9[ 12 . 1 0[) and 
consistents with the previous works [171419] . We also find the ratio of the ion-neutral collision 
frequency to the ion gyrofrequency, a, will directly affect the heating process. When the ratio 
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Figure 2: The temporal evolution of the perpendicular kinetic temperatures normalized with 
respect to their initial values Tj. The various input parameters, (SB W 2 / B 2 , v p /va) = (0.05, 
0.07) and (0.12, 0.07) are shown in Fig.2(a) and Fig. 2(b), respectively. 
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a is larger, the heating process becomes less efficient. However, the low-beta plasmas are 
still significantly heated, as indicated in Fig. 2 where the temperature T±(t) is much larger 
than the initial temperature Tj. This phenomenon may provide a partial explanation as to 
why the temperature of the chromosphere is higher than the photosphere while lower than 
the corona, as the fraction of ionized particles increases from the photosphere through the 
chromosphere and to the corona [27]. Furthermore, we can see from the evolution of the 
population in Fig. 1(e), (f), (g), and (h) that there exists a dynamic equilibrium between 
the heating and cooling of the ions, which is reached after a period of time. This is also 
demonstrated by the oscillations in temperature which decrease over time as shown in Fig. 2 
for the cases a ^ 0. 

The results of the test particle calculations are strongly consistent with our analytical 
predictions. The ratio of frequencies a in the analytic solution Eq. fj2.9j) is in the term e~ Uint 
where v in = af2 ; therefore, a larger value of a will result in a lower value of e~ Vint , which 
will affect the value of ion kinetic temperatures as shown in Fig. 2. Given the parameters 
u{t) ^0, Vp = 0.07v A , 8B 2 W /Bl = C = 0.05, a = 0.1 and t = 40, the solution of Eq.O can 
be reduced to 

| u± | „ Efc [VA - V\\(0)] (Aq + Win) c - ik(vAt - z )- iv)k 

^ y/CK-^O)] (1 + to) ik{vAt _ z) _ i(pk 

1 + a 2 

~ \/C x v A | (1 + ia)e~ tk(vAt - z) -^ k | w x v A VT+a? w 3.2v p (2.12) 

where we use the approximations e~ £ ~ when e > 3, \va — v»(0)\ ~ \va\ and 1 + a 2 ~ 1. 
Thus the stable value of velocity is 3.2v p which is in accordance with the ring structure of 
velocity distribution shown in Fig. 1(h). 

The physical mechanism of this heating can be described as follows: in low-beta plasma, 
nonresonant wave particle scattering by Alfven waves can lead to randomization of the par- 
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ticle motion transversely to the background magnetic field and thus effectively to heating of 
the ions. Due to the fact that the neutrals are given an approximately invariable bulk veloc- 
ity, the heated ions will lose energy from ion-neutral elastic collisions throughout the process. 
This leads to an overall reduced efficiency of the ion heating relative to the simulation with 
no ion-neutral collisions. 

III. CONCLUSION 

In summary, we show that partially ionized low beta plasmas can be heated by a spectrum of 
low-frequency (0.0014 Hz ~ 0.0079 Hz) Alfven waves with large amplitude (2.6 kilometers per 
second in amplitude), which have been observed in the low-beta (where /3 reaches a minimum 
with (3 < 3 x 10 -2 [3]) chromosphere PQ. This is contrary to the linear theory, according 
to which ions can only be heated through resonant interactions with Alfven waves. In our 
model, the frequencies of the Alfven waves are much lower than the ion cyclotron frequency, 
so the cyclotron resonant condition is not met. This also ensures that the waves will not be 
damped by ion- neutral collisions [2"3~H2"o] . We find that the amount of heating depends on the 
ratio a, where less heating occurs when the importance of ion-neutral coupling is increased. 
However, in all cases, significant heating of the ions occurred. Furthermore, we show that the 
velocity distribution will form a ring structure during the heating process when elastic ion- 
neutral collisions are considered, which is consistent with our analytic approximation. Our 
results demonstrate significant heating of ions through nonresonant Alfven wave interactions 
in partially ionized plasmas, which may provide an alternate source of heating in the solar 
chromosphere than previously considered. 
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